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We present a generalized mode-coupling theory for a dense binary fluid mixture. 
The theory is used to calculate molecular-scale renormalizations to the stress-tensor 
autocorrelation function (STAF) and to the long-wavelength zero-frequency shear 
viscosity. As in the case of a dense simple fluid, we find that the STAF appears to de- 
cay as t~ 3 / 2 over an intermediate range of time. The coefficient of this long-time tail 
is more than two orders of magnitude larger than that obtained from conventional 
mode-coupling theory. Our study focuses on the effect of compositional disorder on 
the decay of the STAF in a dense mixture. 



1 Introduction 

In a simple fluid the five conserved densities relax very slowly at small wave vectors. 
The nonlinearities in the hydrodynamic equations provide a further mechanism for 
slowing down the decay of fluctuations and are responsible for the long-time alge- 
braic tails of the correlation functions that determine the transport coefficients [1]. 
The coefficients of the algebraic tails arising from these long- wavelength or conven- 
tional mode-coupling effects are, however, often much smaller than those observed by 
molecular dynamics. In particular, the conventional mode-coupling theory yields a 
coefficient of the long-time tail of the stress-tensor autocorrelation function (STAF) 
of a dense hard-sphere fluid which is 500 times smaller than observed in simula- 
tions [2]. It has been shown that this difference can be accounted for in terms of 
finite- wave- vector or generalized mode-coupling effects [3]. In a dense fluid density 
fluctuations relax slowly even on molecular length scales. As a consequence, density 
nonlinearities lead to a renormalization of the transport coefficients via a finite- wave- 
vector mode-coupling mechanism [4,3] that can qualitatively account for the size of 
the coefficient of the long-time tail of the STAF [2] for the shear-dependent viscosity 
[5,6] observed in computer simulations. In addition they provide a good description 
of the behavior of the velocity autocorrelation function at intermediate times [7]. 
A self-consistent implementation of the generalized mode-coupling theory also pre- 
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diets, at a critical density, the transition to a glassy state where structural relaxation 
is frozen [4,8-10]. On the other hand, a more detailed analysis has shown that the 
transition is an artifact of the approximation used and is cutoff when additional 
couplings are included [11]. The generalized mode-coupling theory has, however, 
provided considerable insight into the dynamical properties of dense simple liquids. 
In particular, it gives an adequate description of the slowing down of the dynamics 
of dense fluids of moderate viscosity, in a regime above the glass transition were the 
relaxa- tion is highly cooperative and nonexponential [12,13]. 

In this paper we are interested in the mode-coupling contribution to the shear vis- 
cosity and in the long-time behavior of the stress-tensor autocorrelation function 
of a dense binary fluid mixture. In a recent paper [14], Erpen- beck reported 
molecular-dynamic studies of an isotopic hard-sphere mixture at moderate densi- 
ties and concluded that in this case the size of the coefficient of the long-time tail 
of the time-correlation function for the shear viscosity observed in the simulations 
agrees with the prediction of conventional mode-coupling theory. In mixtures, as 
in one-component fluids, the conventional mode-coupling contribution to the STAF 
only accounts for the kinetic part of the correlation function. The kinetic contribu- 
tion dominates at low and moderate densities and in this regime the conventional 
mode-coupling theory agrees with the simulations. Even in one-component fluids 
[15] it is only at high densities (na 3 > 0.7) that the largest contribution to the 
correlation function is the potential one and the observed coefficient of the t~ 3 ^ 2 tail 
of the STAF is much larger than predicted by conventional mode-coupling theory. 
The situation should not be different for mixtures. The agreement of the simulations 
with conventional mode-coupling theory can only occur at the relatively moderate 
densities considered by Erpenbeck. His observations do not preclude the possibility 
of obtaining large finite-wave-vector mode-coupling corrections to the stress-tensor 
autocorrelation function at higher densities. To our knowledge such effects have not 
yet been investigated in computer simulations of mixtures. 

Here we evaluate the generalized (to include finite- wave- vector effects) mode-coupling 
contribution to the STAF and to the shear viscosity of a binary mixture of hard 
spheres of different sizes and masses. The generalized mode-coupling theory takes 
into account that in a dense viscous liquid structural rearrangements are very slow 
due to a close packing of the molecules. The slow dynamics of density fluctuations 
on molecular length scales affects the macroscopic transport properties through a 
nonlinear mode-coupling mechanism [4,3]. At large wave vectors the most important 
of these contributions involves the product of two density fluctuations [16,3]: this 
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contribution is the slowest decaying and has the largest amplitude. In the case of 
a one-component fluid, the dynamics slows down significantly at only one value of 
the wave vector corresponding to the location of the first maximum of S(k) [17]. 
In contrast, in a dense mixture of spheres of different sizes fluctuations in the mass 
densities of both species can be long-lived on molecular length scales [18] and the 
dynamics can slow down over a range of wave vectors. In paper I [18] we obtained 
a set of linear generalized hydrodynamic equations that describe the dynamics of 
the density fluctuations of the two species down to molecular scales. We found that 
when the packing fractions of the two species are comparable and all three partial 
static structure factors of the mixture are peaked, the density fluctuations of the two 
species decay very slowly at two length scales corresponding to the molecular diam- 
eters of the two species. We now consider the nonlinear generalized hydrodynamic 
equations and discuss the renormalization of the transport coefficients by density 
nonlinearities. The interplay of the density fluctuations of the two components leads 
to finite-wave-vector mode-coupling corrections to the correlation functions which 
determine the transport coefficients. 

We are interested in a binary mixture of hard spheres as the simplest model where 
one can study in some detail the role of compositional disorder in slowing down the 
fluid dynamics. Compositional disorder seems to be an essential feature of the sim- 
plest glass formers (metallic glasses). Our work is focused towards understanding 
how the structural relaxation of the two species and interdiffusion affect the dense 
fluid dynamics. We choose to study the mode-coupling correction to the stress- 
tensor autocorrelation function which determines the shear viscosity as the simplest 
of the various mode-coupling effects in a mixture. 

The self-consistent implementation of the generalized mode-coupling theory of bi- 
nary mixtures has been discussed recently by several authors as a model of the glass 
transition [19,20]. Our work differs from this work on mixtures because here we 
simply study the mode-coupling corrections to lowest order in perturbation theory, 
without attempting any self-consistent closure. We do so for two reasons. First we 
believe that one needs to understand well the simplest perturbation theory before 
making it self consistent. Second, and most importantly, it has been shown [3] that 
for one-component Ouids a naive mode-coupling theory of the type we discuss here 
describes well the slowing down of the dynamical properties above the glass tran- 
sition, in a regime of moderate viscosity where the dynamics is cooperative (and 
therefore nonexponential), but not yet activated [12]. Our objective is to put the 
description of the dynamics of dense mixtures in this regime on the same solid foot- 
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ing as that of one- component fluids. 



The outline of the paper is as follows. In Sec. II the. basic theory used here is 
described. The perturbation theory in the nonlincarities in the hydrodynamic equa- 
tions is then outlined in Sec. Ill, where an expression for the renormalized shear 
viscosity is obtained. In Sec. IV we present our results for the STAF and the shear 
viscosity for hard-sphere fluid mixtures. We analyze the concentration dependence 
of the STAF for a mixture and conclude the paper by discussing our results and by 
making some remarks about connections to earlier theoretical studies and computer 
simulations. 



2 Nonlinear generalized hydrodynamics 

The hydrodynamic properties of a binary fluid mixture are described in terms of 
the six conserved densities: the mass densities pi(r,t) and p 2 ( r ?^) P-( r )t) of the two 
species, the total momentum density g(r, t), and the energy density e(r,t) [18]. We 
begin by writing a set of hydrodynamic equations that have the same structure both 
at large and short (molecular) length scales. The starting point for the derivation 
of these generalized hydrodynamic equations is the free-energy functional for an 
inhomogeneous equi- librium binary fluid mixture, given by 

dr\^ + f]+F v . (2.1) 
_2pi 2p 2 _ 

The potential part Fy of the free-energy functional can be written exactly in the 
form of an expansion in the fluctuations Sn s (r, t) —no s of the number density n s (r, t) 
of species s, with s = 1, 2, from its equilibrium value n 0s , [21] 

5F V k B T 



\n{\ 3 Ts n s (r))-Yl I dv'C ss ,{\r - v'\)8n s ,{r' ,t) 

s'=l J 

I r ' I dr"C ss t s »(r,r',r")6n s r(r',t)6n s »(r",t) + 
s '=i/=r •* 



Sp s m 

2 2 



• (2-2) 



Here m s is the mass of a particle of type s, T is the temperature, Xt s = (h 2 /2irm s kBT) 1 / 2 
is the thermal wave-lngth of particles of type s and C ss '(|r — r'|) and C ss r s rr(r, r'r") 
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are the equilibrium two-particle and three-particle direct correlation functions. From 
this free energy, nonlinear fluctuating hydrodynamic equations can be derived by 
the standard Poisson bracket method outlined in Appendix A [22]. Here we neglect 
temperature fluctuations since they do not play an important role in determining 
the renormalization to be discussed in Sec. III. The equations for the partial mass 
densities of the two species and the total momentum density are 



<9 t pi(r, t) + nabla 
d t p2{r,t) + nabla 



P 



P2 
P 



g(r,t) 



-G x {T,t) + h{v,t), 
-G 2 (r,t) + f 2 (r,t), 



(2.3) 
(2.4) 



d t g a (r,t) + 5 a 



9a9f3 



+ ^2 p s d a 



6F, 



V 



8p s 



1=1 

-G ga (v,t) + f ga (v,t). 



(2.5) 



Here Gi(r, t) and G2(r, t) account for diffusive dissipation of species type 1 and 2, 
respectively, and G 9a (r, t) corresponds to viscous dissipation. We write these terms 
in the form 



G s (v,t) = d a J2 J ^7S(r,r'; 

s'=l 



n s )& 



5F 



SpAr') 



for s — 1, 2, and 



G ga (r,t)= dv'Y 



5F 



a(3\ 



(2.6) 



(2.7) 



From the equation of continuity for the total density fluctuations p(r, t) = pi(r, t) + 
P2(v,t), it follows that Gi(r, t) + G 2 (r?£) = 0. As we have seen in /, the only fluc- 
tuations that are slowly varying on molecular length scales are the fluctuations of 
the partial densities of the two species. Large-wave-vector momentum fluctuations 
relax on short-time scales and momentum behaves like a nonhydro dynamic (noncon- 
served) variable. We have therefore written Gi(r, t) and G^ 1 "; t) in a form explicitly 
consistent with the general structure of dissipative terms pertaining to conserved 
variables. The form of G 9a (r, t) corresponds to that of dissipation associated with 
a non-conserved variable [23]. On the other hand, it should be kept in mind that 
our dissipative kernels %s* r Q/ 3 are nonlocal and reproduce the well-known transport 
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coefficients in the long- wavelength limit. 

In the above equations /i(r, t), /2(r, t), and f 9a (r,t) are Gaussian random fluc- 
tuating forces which are related to the corresponding linearized dissipative kernels 
through the fluctuation-dissipation theorem. The random forces have zero average 
and correlations given by 

(f s (r,t)fAr\t')) = 2k B T 1 ^(r,r')d a d^(r-r')6(t-t'), (2.8) 

(f 9a (v,t)f gg (r\t')) = 2k B TT aP (r,r')5(v - r')5(t - t>). (2.9) 

Equations (2.3)-(2.5) follow from general considerations if one neglects temperature 
fluctuations and non-Markovian effects. 

To complete the equations, one needs to specify the dissipative kernels. We will 
neglect dissipative nonlinearities in our analysis. For a dense one-component fluid 
a theory that only retains thermodynamic nonlinearities predicts a t -3 / 2 tail for the 
STAF with a coefficient that is a factor of 2 smaller than that observed in simu- 
lations [7]. The inclusion of dissipative nonlinearities accounts for this factor-of-2 
difference [17]. For simplicity here we neglect dissipative nonlinearities and assume 
r Q/3 (r, r';n s ) ~ T a/3 (r, r'; n s0 ) and j^(r,r';n s ) ~ j^, (r , r' ; n s0 ) . The simple fluid 
results indicate that even with this approximation the theory will provide a semi- 
quantitative description of the behavior of the STAF at intermediate times. We also 
neglect convective non-linearities which are responsible for the asymptotic t~ 3//2 tails 
in the STAF. At high density the coefficient of the asymptotic tail is much smaller 
that the contribution discussed here. 

It is convenient for the discussion below to consider the Fourier transform of the gen- 
eralized hydrodynamic equations (2.3)-(2.5). The Fourier components of momentum 
fluctuations are separated into longitudinal and transverse parts, 

g(k,t) = k#(k,t) + K_gnQL,t), (2.10) 

i=l,2 

where k = k/fc is a unit vector along k and k^ are unit vectors normal to k. After 
neglecting convective and dissipative nonlinearities and retaining only quadratic non- 
linearities in the density fluctuations, the Fourier-transformed equations are given 
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by 



d t p s (k,t) + ik P ^-g l (k,t) = V" 7,./ ~ tvCW'W] ^"(k' ^ + f s (k,t), 

(2.11) 

for s = 1, 2, and 

2 

[<% + 7i(fc)]^(k,t) + zfc ^ i S ss> ~ n 8 C a8 >(k)]5 nal (k,t) 

s,s'=l 

= I L (k,t) + F L (k,t) (2.12) 

[dt + Mk)\9T i (Kt) = / Ti (k,t) + / r< (k,t). (2.13) 
In writing the above equations we have used that for an isotropic fluid 

r a/3 (k) = p [k a k/3^L(k) + (6 a/ 3 - k a kp)y T (k)], (2.14) 

% s ,(k) = 5 af37ss ,(k). (2.15) 

In the long- wavelength-limit the viscous kernel r a/ g(k) reduces to the usual combi- 
nation of gradients and bare shear and bulk viscosities rj s and £, 

limr Q/3 (k) = fc a A^(C + ~r) a ) + (& 2 £ Q/ 3 - k a kp)r] s . (2.16) 

The diffusion kernels j S s'(k) are related to the diffusion constants D ss > of a binary 
mixture, according to 

£> ss ' = lim 5 - r-y(fr)- (2-17) 
fc^o y/m s p s m s <p s < 

Finally, the nonlinearities in Eqs. (2.12) and (2.13) are given by 

2 

4M = |E / 7Si[^(k,q) + K'(q ) k)] 

s,s'=l 

x5 ns (q,t)cL,(k-q,t), (2.18) 
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x8 ns (q,t)8 n ,(k-q,t), (2.19) 



with 



V2(k,q) + K?(q>k) = ^[C ss '(|k - q|) - C. y (g)], (2.20) 



V£(k,q) + y£(q,k) = ^[C ss ,(|k-q|)-C ss ,( ? )], 

+ ^^^"^"(q-kj-q). (2.21) 

P s"=l 

The linearized form of Eqs. (2.3)-(2.5) is identical to the linear equations obtained 
in paper I for a binary fluid of hard spheres from the Enskog theory. For the case of 
hard-sphere fluids one can derive explicit analytical expressions for the dissipative 
kernels. These expressions were obtained in I. The viscous kernels are given by 



lL\k) = — > -f- di(fcow), 2.22 



7r(fc) = 2^ T d 2( k<J ss'), (2.23) 



3p ^-f *B , 

where £e ss ,, is the Enskog mean free time between collisions, 

1 Ay/7rn s n s >ag S , 

= (2/W) 1 / 2 %SS '' (2 ' 24) 

with /i ss / = m s m' s /(m s + m s i) and c/ii ss / the pair correlation function of species s and 
s' evaluated at contact, chi ss / = g ss >{r = cr ss r). The functions d\(x) and d 2 (x) are 
defined in terms of the spherical Bessel functions ji(x) as d\(x) = 1 — jo(x) + 2j 2 (x) 
and o^O^) = 1 — jo{x) — 32{x). The four ditTusion kernels 7 ss '(fc) are related via the 
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Onsager reciprocity relations, which require ■yn(k) = 722 (k) = —'yilik) = — 721 (fc). 
Therefore there is only one independent diffusion constant in a binary mixture. For 
hard-sphere fluids it is given by 

7ii (k) = 2 

where 

d s Ax) = ^ - (-l) s+s '[j (x) - 2j 2 (x)} (2.26) 

Ps' 

This concludes the derivation of the nonlinear generalized hydrodynamic equations. 
In the next section we describe the perturbation theory to evaluate the mode- 
coupling contribution to the STAF. 



2fi a 



■f y/m s m sl p s p s ftE 3sl 



d ss '(ka s 



(2.25) 



3 Mode- coupling theory 

We are interested in evaluating the generalized mode-coupling contribution to the 
coefficient of shear viscosity r] s and to the corresponding STAF, p v (t). They are 
related by the Green-Kubo formula, 

POO 

Vs=(3 dtpr,(t). (3.1) 
Jo 

The shear viscosity governs the decay of the transverse momentum fluctuations, 
Cr(k, t) = (<7T(k, t)gT(— k, 0)). When nonlinearities are neglected, the generalized 
hydrodynamic equations discussed in Sec. II give 

[d t + lT (k)}CP(k,t)=0, (3.2) 

where 

lim7 T (A;) = k 2 r] E /p , 

with t]e the Enskog shear viscosity of a mixture [24,25]. The solution to (3.2) is 
obtained immediately, 
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CP(k,t) = e-^ (k)t C T (k,0) 

~e VEk2t/p0 C T (k,Q), (3.3) 

where the second equality applies in the limit k — > and CV(/c,0) = p/fi. When 
quadratic nonlinearities are retained in Eq. (2.13), we obtain 



[d t+jT (k)]C T (k : t) = ~ ]T | / ^[K^(k,q)+Kl(q,k)](5n s (q,t)^(k-q,t)^(-k,0)), 

(3.4) 

where the prime on the integral sign denotes a short distance cutoff. The effect of the 
nonlinearities on the right-hand side of Eq. (3.4) can be approximately incorporated 
by a perturbation theory that reexpresses the correlation function on the right-hand 
side of Eq. (3.4) in terms of CV(k, t). The theory is described in Appendix B. 
Equation (3.4) is then recast in a form where the contribution from nonlinearities 
enters as a non-Markovian renormalization of the transverse damping 7^, 



[d t + J T (k)]C T (k, t) = f dr VVk, t - r)C T (k, r) = 0, 

^0 



(3.5) 



with 

E(M) = ^E f ^\^M\ 2 eM^^]t (3.6) 

and 

2 2 

|^(k,q)| 2 E [Kl'(k,q)+V£(q,k)K s ^^^ 

s,s'=l s",s"'=l 

(3.7) 

Explicit expressions for the amplitudes a^ s , are given in Appendix B. 
It is convenient to consider the Laplace transform of (3.5), 

[z + <y*{h, z)]C T (k, z) = C T (k, t = 0), (3.8) 



where 



poo 

C T (k,z)= dt e- zt C T (k,t) (3.9) 
Jo 
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for Re(z) > 0. The renormalized transverse dissipative kernel 7^ is the sum of its 
bare part 7t(&) and a mode-coupling contribution ^7™° 

^(K,z)= lT (k) + 5^ c (k,z), (3.10) 

with 

5^ c (K,z) = dte- 2l V(k,i). (3.11) 
Jo n 

The long-wavelength limit of the transverse dissipative kernal 7^ (k, z) defines the 
frequency-dependent shear viscosity r)f(z) according to 

hin^{k, z) = k 2 ri*{z). (3.12) 

The mode-coupling contribution to the renormalized viscosity, r)f"(z) = 7]E + 5ri™ c (z), 
is then given by 

^rW = lim-/ dte-*£;(M). (3.13) 
jo ^ 

By comparing Eqs. (3.1) and (3.13), we see that the long-wavelength limit of 
5^_(k, t) also determines the mode-coupling contribution to the STAF p v (i) ac- 
cording to 



^ cn (t) = liml^(k,t). (3.14) 

v 

Finally, 5p° cn (t) and 5rjf c are explicitly given by 



W = \ E r?03lV(O,q)| 2 

Xe fe(q)+^(-q)]t ? (3. 15 ) 

and 

^ r --E r iV^'\r (3-16) 

/s 2 ^7 (2tt)3 [^(«j)+^(-q)] 
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4 Results 



In this paper we have used a generalized mode-coupling theory that takes into ac- 
count molecular scale effects to compute the STAF of a dense binary-fluid mix- ture. 
As in the case of a one-component fluid of packing fraction 77 = (ir/6)na 3 = 0.46, 
in a dense mixture the coefficient of the long-time tail of the STAF evaluated from 
generalized mode-coupling theory is about two orders of magnitude larger than the 
corresponding coefficient obtained from conventional mode-coupling theory. 

Here, we consider a mixture of hard spheres of sizes (J\ and er 2 , with er 2 > <j\. As 
discussed in paper I, the slowing down of the dynamics of the dense mixture on short 
length scales is directly related to the short-range correlations that characterize the 
fluid structure and that show up as peaks in the partial static structure factors. 
These are in turn determined by the relative packing fractions of the two species, 
rji = im\o\/§ = r]Xia 3 /(x 2 + Xi« 3 ) and r] 2 = 7rn 2 cr|/6 = 77£ 2 /(x 2 + Xi« 3 ), where 
77 = rji + r] 2 is the total packing fraction and xi = n\/n and x 2 = n 2 /n, with 
n = Hi + n 2 , are the concentrations. Notice that in our calculation we keep 77 
constant 77 = 0.46 in all the data discussed below), while changing the concentration 
x 2 of the large spheres. Finally, in this paper we only present results for mixtures of 
spheres of equal masses. This is done mainly to dissociate the role of the structural 
parameters from the role of the mass ratio in slowing down the dynamics. All times 
are normalized to an Enskog time given by 



1 4^(77,1 + n 2 )cr 



t E - (2M» * 12 ' (41) 
This time scale is neither the Enskog mean free time of a hard-sphere fluid of type-1 
nor type-2 spheres. It is chosen because it depends only weakly on the concentration 
of each species. 

We first discuss the conventional mode-coupling contribution to the STAF of a mix- 
ture, denoted by 5p™ nv (t). As in one-component fluids, one finds Sp^ onv (t) ~ At" 3 1 2 
[26]. This algebraic tail arises from the slow decay of fluctuations in the conserved 
densities at long wave-lengths and the feedback from these long-lived fluctuations 
through the convective nonlinearities in the hydrodynamic equations. It describes 
the true asymptotic long-time behavior of the correlation function. The coefficient 
A was evaluated many years ago by Pomeau [26] and, as in one-component fluids, 
it was found to be very small. In Fig. 1 we show 5p™ nv (t) = / '77^ / 've)^ r {t) 
as a function of (t/tE)~ 3 ^ 2 for three values of the concentration x 2 = 77, 2 /n of larger 
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particles. Here t]e is the Enskog shear viscosity of the mixture [24]. The coefficient 
of the long-time tail, given by the slope of the straight lines of Fig. 1, has a weak 
dependence on the concentration x^- The slope is also quite independent on the 
total packing fraction r\. 

The conventional mode-coupling contribution to the STAF should be compared to 
the generalized mode- coupling contribution in Eq. (3.15) and denoted here by 
5p^ cn (t). Equation (3.15) can be written as a one-dimensional integral over the 
magnitude q of wave vector. This integral has been evaluated numerically using an 
upper cutoff g m axO"i2 = 20. We have checked that the contribution to the integral 
from reduced wave vectors larger than 20 is negligibly small. In Fig. 2 we show 
the generalized mode-coupling contribution Sp^ en (t) as a function of (t/t E 3 ^ 2 for a 
packing fraction 77 = 0.46. The curves are well approximated by straight lines over a 
range of time scales, indicating that the generalized mode-coupling contribution to 
the STAF appears to decay as ~ t~ 3 ^ 2 in this region. The generalized mode-coupling 
contribution evaluated here describes the decay of the correlation function over the 
intermediate time scales that are generally accessible to computer simulations. It is 
very hard to probe by simulations the truly; asymptotic long-time behavior of the 
correlation functions described by the conventional long-time tails. 

The apparent failure of our curves in Figs. 1 and 2 to reduce to the result obtained 
for a one-component fluid in the limit of either small (X2 = 0.01) or large (X2 = 0.9) 
concentration of large spheres is only an artifact of our choice of the units of time 
[here we have scaled time with given by Eq. (4.1)]. If in each case (X2 << 1 and 
X2 — 1) we scale the time with the Enskog mean free time of the majority compo- 
nent, tEii and tE22, respectively [see Eq. (2.24)], then our results for X2 = 0.01 and 
X2 = 0.9 are essentially identical to those for a one-component fluid. 

By comparing Figs. 1 and 2 we see immediately that at the large values of the 
density considered here the generalized mode-coupling contribution to the STAF is 
about two orders of magnitude larger than the conventional mode-coupling contri- 
bution (the vertical scales in Figs. 1 and 2 differ by a factor of 10 3 ). We also find 
that Sp^ en (t) depends strongly on the fluid density. For 77 = 0.23 (the value used in 
Erpenbeck's simulations) 5p^ cn (t) is smaller than 5p™ nv (t). This is consistent with 
the fact that the generalized mode-coupling contribution accounts for the potential 
part of the correlation function that dominates at high density. One of the factors 
that controls its size is the height of the first peaks of the static structure factors of 
the fluid, which in turn are large only at high density. 
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The generalized mode-coupling contribution given in Eq. (3.15) is the sum of four 
terms, corresponding to the coupling of two density modes (- -), the coupling of 
two diffusion modes (+ +), and two equal terms corresponding to the coupling of 
one diffusion mode and one density mode (+ - and - +). In Figs. 3-5, we show 
separately these contributions to the STAF as well as their sum for three values of 
the concentration x 2 For a small (#2 = 0.01, Fig. 3) and a large (x2 = 0.9, fig. 
5) concentration of large spheres the generalized mode-coupling contribution to the 
STAF is dominated by only one of these four terms. For X2 << 1 the density mode 

(z—) describes the relaxation of total density fluctuations which are mainly deter- 
mined by the majority component (type 1, in this case). The coupling of two of 
these density modes gives then the dominant generalized mode-coupling contribu- 
tion to the STAF. For x% ~ 1 it is the mode we label diffusion mode (z+) that 
describes the relaxation of density fluctuations and the coupling of two of these 
modes dominates the mode-coupling contribution to the STAF. At an intermedi- 
ate value of the concentration of large spheres (x% = 0.2, Fig. 4) we find that the 
contributions from all four terms in Eq. (3.15) are of comparable magnitude. For 
large or small values of the concentration of one of the two species (here the large 
spheres), the mixture resembles a one-component fluid and the density fluctuations 
of the majority component govern the slow structural relaxation of the fluid. The 
effect of diffusion of the low-concentration component on slowing down the decay of 
the STAF appears to be small. For the intermediate value of the concentration the 
packing fractions of the two species are comparable and the mixture structurally 
resembles a metallic glass. The mixture is structurally very disordered and all three 
partial static structure factors are peaked at well-separated values of the wave vec- 
tor, determined by the sizes of the two species. As discussed in paper I, fluctuations 
in the density of both species decay slowly at molecular length scales. The two 
extended modes of diffusion and total density fluctuations cannot be decoupled: 
both soften appreciably at large wave vectors and contribute to the mode-coupling 
integral. There is therefore a range of large wave vectors where density fluctua- 
tions decay slowly. This should be contrasted to what happens in the case of a 
simple fluid, where the mode-coupling integrals are dominated by the single value 
of the wave vector corresponding to the first maximum of the static structure factor. 
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We now discuss in some detail the effect of compositional disorder on the decay of 
the STAF in a dense mixture. For the present discussion we normalize all times to 
t* given by 

1 _ Vgniof + n 2 a\) , 2 

^ (2f3fj, 12 y/ 2 ^i^ 11 + X 2X22 + ^1^2X12; • 14. 

This time scale, unlike t E , has the feature of going to the one-component fluid limit 
for x 2 << 1 and x 2 — 1- 

We first consider the case of conventional mode-coupling. In Fig. 6(a) we show 
5p r v mconv *{t) = (f3t*/7] E )5p r v mconv (t) as a function of (t/t*) for three small values of 
the concentration x 2 = n 2 /n oi larger particles. Increasing the concentration of 
large particles in a fluid of small spheres slightly shifts the characteristic time scale 
governing the decay of the STAF to larger times. For a given value of {t/t*), the 
STAF of the mixture is then larger than that of a one-component fluid. This is, 
however, a very small effect and it is entirely due to the concentration dependence 
of the Enskog shear viscosity of the fluid. A similar behavior is observed when small 
particles are added to a fluid of large particles [see Fig. 6(b)]. 

The concentration dependence of a generalized mode-coupling contribution to the 
STAF is quite different. In Figs. 7(a) and 7(b) we display 5pf n *(t) = (f3t* /7] E )5pf n (t) 
as a function of {t/t*) for various values of the concentration x 2 = n 2 /n of larger par- 
ticles. For values of x 2 close to 1 [Fig. 7(b)], when the large spheres are much more 
abundant than the small spheres, increasing the fraction of small spheres results in a 
shift of the characteristic time regime where the mode-coupling contribution to the 
STAF is large to shorter times. The shift is in the opposite direction as compared 
to that observed for the conventional contribution and it is more pronounced. A 
similar behavior is also observed for small values of x 2 [Fig. 7(a)]. 

A comparison of Figs. 6 and 7 shows that the concen- tration dependence of §p& en (t) 
is stronger than that of 5p™ nv (t). This indicates that the degree of compositional 
disorder in the fluid plays a role in the dense fluid dynamics. 

In order to assess more precisely when compositional disorder enhances large-wave- 
vector mode-coupling effects, we have also evaluated the mode-coupling contribution 
to the shear viscosity, Sp^}. In Fig. 8 we show Spf^/r] E as a function of x 2 for a 
few values of the size ratio a. Here 77^ is the Enskog shear viscosity of the mix- 
ture. In order to properly interpret the results displayed in Fig. 8 it is important 
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to first clarify the concentration dependence of the Enskog shear viscosity. At the 
high densities considered here the Enskog shear viscosity of a hard-sphere fluid is 
dominated by the potential contribution. For a one-component fluid this is given 
by Ve = ^\/vrm//5[(r2(T 3 ) 2 x/cr 2 ]. If the reduced density na 3 is kept fixed, then 
~ I/ " 2 and the viscosity decreases as the sphere diameter increases (the pair 
correlation at contact \ only depends on na 3 ). The Enskog shear viscosity tje of a 
binary mixture depends on the concentration x 2 of larger particles and reduces to 
the Enskog viscosity of a fluid of hard spheres of size a±, t]ei, for x 2 = 0, and to 
the Enskog viscosity of a fluid of hard spheres of size a 2 > ai, r]E2, for x 2 = 1, with 
VE2 < Vei- Compositional disorder always increases the bare viscosity in the sense 
that for x 2 « 1, r\ E > r] E1 and for x 2 — 1, t]e > r] E2 . 

In contrast we see from Fig. 8 that the concentration dependence of the generalized 
mode-coupling contribution 5ry™° (relative to t]e) is quite different. The dashed line 
in the figure indicates the value of 5rf^ c /riE for a one-component hard-sphere fluid 
with na 3 = 0.88 as obtained before by other researchers [27]. Our results show that 
in a mixture with a small fraction of large spheres (a colloidal suspension) the gen- 
eralizai mode-coupling corrections are larger than in a one-component fluid. This 
corresponds to the fact that the motion of large spheres in a dense background fluid 
is very slow. Shear relaxation also gets slower as the size difference increases (smaller 
a). On the other hand, in a mixture containing a small fraction of small spheres in a 
fluid of large spheres (x 2 ~ 1) the generalized mode-coupling correction to the shear 
viscosity, while still larger than the one-component value, is considerably smaller 
than obtained in a colloidal-fluid-like mixture (x 2 << 1). This is because small 
spheres can move rather easily through a matrix of large spheres. These considera- 
tions are strictly meaningful only for a = 0.5 and a = 0.7, that is, for mixtures of 
spheres of rather disparate sizes. The a = 0.9 value shows how the curves converge 
to the one-component fluid limit. Our findings are in agreement with the results of 
simulations of mixtures by Mountain and Thirumalai [28]. These authors reported 
that a mixture with a small fraction of small particles can crystallize more easily 
than a mixture with a few large particles, due to the fact that small particles can 
move more easily in a matrix of large particles compared to large particles in a ma- 
trix of small particles. 
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5 Discussion 



In this paper we have formulated a generalized mode-coupling theory for a dense 
binary fluid mixture that takes into account the slowing down of the fluid dynamics 
at molecular scales due to the closely packed fluid structure. This has been used 
to evaluate the renormalization of the STAF due to density nonlinearities in the 
hydrodynamic free energy. As in the case of a one-component fluid [3], in a dense 
fluid mixture the long-time tails in the STAF are about two orders of magnitude 
larger than the conventional long- wavelength contributions. The concentration de- 
pendence of 5p& cn (t) is stronger than that of 5p c ° m (t). This clearly indicates that 
at high densities the degree of compositional disorder in the fluid affects the static 
fluid structure and plays an important role in slowing down the fluid dynamics. In 
this section we make a few comments on the results obtained in this paper. 

(1) In our analysis we have considered a dense binary mixture of different sizes and 
equal masses at different values of the concentration. The packing fraction of 
each species is the main parameter controlling the static structural properties 
of dense fluids at molecular length scales. For this reason concentration and 
size ratio have qualitatively the same effect on the short-wavelength dynam- 
ics. We find that at intermediate values of the concentration when the fluid 
is structurally disordered, there is a finite range of wave vectors over which 
the dynamics slows down and consequently the mode-coupling integrands are 
large over a range of length scales. This is reflected in the fact that all possible 
couplings of the density fluctuations of the two species make comparable con- 
tributions to the total 5p^ cn (t). The presence of a range of length scales and 
hence a distribution of time scales for a disordered dense binary mixture that 
we predict is a characteristic feature of glass-forming liquids [29]. In contrast, 
in a simple fluid the mode-coupling integrands are very sharply peaked at the 
position of the first maximum of the static structure factor. The presence 
of compositional disorder shifts the time regime where the generalized mode- 
coupling contribution to the STAF is large to shorter times. In order to test 
our results one needs to perform molecular-dynamics simulations for mixtures 
of dissimilar sizes at high densities. 

(2) In this paper we only considered mixtures of spheres of equal masses. This 
choice was motivated by our desire to focus on the role of frustration caused 
by size difference in slowing down the fluid dynamics. In addition, both our 
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previous work [18] and computer simulations [20] have indicated that the slow- 
ing down of the dynamics in a dense, closely packed binary mixture depends 
more weakly on the mass ratio than on the size ratio of the two species. 

(3) The nonlinear fluctuating hydrodynamic equations presented here can also 
be used to discuss the mode-coupling theory of the glass transition in fluid 
mixtures. In this case the mode-coupling contribution to the time- depen- 
dent density-density correlation function is evaluated self-consistently. The 
self-consistent mode-coupling theory of mixtures has been discussed by other 
authors [19]. It involves a considerable amount of computation that tends to 
obscure the physics. It would be interesting to consider a simplified model 
obtained, for instance, by replacing the static partial factors of the mixture 
by constant values with appropriate wave-vector cutoffs or by Lorentzians. A 
calculation of this type has been carried out by Mazenko and Vails [30] for a 
one-component fluid. 
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Appendix A 

In deriving the nonlinear tluctuating generalized hydrodynamic equations of a bi- 
nary tluid mixture we follow the approach described by Ma and Mazenko [22] . The 
Langevin equation governing the time evolution of the hydrodynamic variables de- 
noted by ifii may be written formally as 



5Mt) = V&] - /.,,—— + Mt) . (Al) 



5F 

'<%(*) 

The reversible part of the dynamics is given by 



W1 = ««Mtt£t-sM- (A2) 

<%(t) (3 5^ jit) 

where Qy[^] = is the Poisson bracket between the variables ipi and F[ip] 

is the free-energy functional given in Sec. II. The damping coefficients Ly and the 
thermal noise fi pertain to the dissipative part of the dynamics. The noise fi obeys 
the fluctuation dissipation relation 

(fi(t)fj(t')) = 2k B TLij6(t - . (A3) 

For our binary mixture the relevant hydrodynamic variables are {ipi} = (p s , g). The 
Poisson brackets needed in the hydrodynamic equations are 

{p s (r),^(r')} = -dj{v - i> s (r') , (A4) 
{g a (r,g a (r')} = d'J(r - r')^(r) - fy(r - r)^(r') . (A5) 
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Appendix B 



In this appendix we brietly describe the perturbation theory that leads to Eq. (3.5). 
In order to close Eq. (3.4) we need to construct an equation for the correlation 
function of three phase functions appearing on the right-hand side of Eq. (3.4). 
This is done most simply by using the invariance of equilibrium correlation functions 
under time translations, 

(8n s (q, t)8n s >(k - q, t)g T (-k, 0)) 
= (8n s (q, 0)5n s ,(k - q, 0)<7 T (-k, -t)) . (Bl) 

By letting t — > — t and k — > — k in Eq. (2.13), multiplying the resulting equation 
with two density fluctuations, and taking the ensemble average, we then obtain 



[~d t + j T (k)](8n s (q, 0)5n s ,(k - q, 0)g T (-k, -t)) 




(^(q , ,-t)^(-k-q,-t)5n s (q,0)5ny(k-q,0)) . (B2) 

We then perform a cluster expansion of the four-point correlation function on the 
right-hand side of Eq.(B2) and neglet cumulants. With this approximation Eq.(B2) 
becomes 

[d t + j T {k)]{Sn s (q, 0)<W(k - q, 0)<? T (-k, -t)) 
A E /^[^(-k.qO + ^v^q'^k)] 

s",s"'=l J ^ ' 

(8n s n(cf,t)8n s (q, 0))(8n a m(-k- q, t)8n s >(k- q, 0)) . (B3) 

The time dependence of the density- density correlation on the right-hand side of 
Eq.(B3) is now evaluated from Eqs.(2.11) and (2.12) by neglected nonlinearities. 
The most important contribution to the mode-coupling integrals comes from wave 
vectors near the peaks of the static structure factors. Momentum fluctuations of 
such large wave vectors decay rapidly and one can neglect the time derivative in 
Eq.(2.12) as compared to the viscous damping 7^ in treating all except the largest 
time scales. The longitudinal momentum can then be eliminated between Eqs.(2.11) 
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and (2.12). The resulting two coupled linear equations for density fluctuations are 
given by 



d t Sn s (q,t) + q 2 b ss ,(q)5n s ,(q,t) = 
where summation over s' is intended, and 

, , \ n s g s ,(q) 7 SS , f 8 ,(q) 
Oss'{q) = — — + 



with 



9i(q) = ^yt 5 ' 22 ^) ~~ V n 2/niS 12 (q)] 

92(g) = ^yt^n^) - V n i/ n 2Si2{q)} 



(B4) 



(B5) 



(B6) 



A(g) 
1 



S 22 {q) + —V n i/ n 2Si2{q) 

m 2 

Sn(q) + —\Jn 2 /n 1 S 1 2(q) 
mi 



with A(g) = Sn(q)S 22 (q) — Sf 2 (q). The approximate modes obtained from Eqs.(B4) 
were discussed in paper I, where it was shown that this simple approximation is 
indeed accurate at large wave vectors. The time decay of the density-density corre- 
lation function is then governed by two relaxation rates that we label density mode 
[z — (q)] and diffusion mode [z + (q)], 



2 2 
z±{q) = |-(6n + b 22 ) ± |-[(& n - b 22 f + Ab^bn] 1 ' 2 . 

The density- density correlation functions are given by 

(Mq,t)<W(-q,0))= J2 <A^y M \ 



(B7) 



(B8) 



where 



a 



Is' 



«2 S '(q) = v 



v /riira s /5'i g /(g)(&22 + z^) - ^n 2 n s >S 2s >(q)bi 2 

[(6n- 6 2 2) 2 + 46 12 6 2 i] 1 /2 
^n 2 n s ,S 2s >(q)(b u + z^) - y/nin s >S ls >(q)b 2 i 
[(6ii- 622) 2 + 46 12 6 21 ] 1 /a 
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(B9) 



The partial static structure factors are defined as {Sn s (q)Sn s i(—q}) 

= ^n s n s 'S ss '(q). Finally the eigenvalues given in Eq.(B7) can be rewritten in the 

more transparent form presented in paper I and rewritten here for completeness, 



z±{q) = -\ 



, 7n(g) J7i(g) | /2(g)] 

\m 1 pi m 2 p2 J . 
1 2 



±- 



1 7n(g) 

pxAghdq) 13 



mipi m 2 p2 



47n(g) 
pf3 2 lL(q) 



fi(q) /2(g) 



1 2 



mi 



m 2 



1/2 



(BIO) 



Here Xt(q) is the generalized isothermal compressibility, defined as 

Sii(q)S 22 (q) ~ Sl 2 (q) 



Xr(q) 



X°T- 



x 2 S n (q) + XiS 22 (q) - 2 v fr[x 2 ~S l2 (q) 
with chi\ — (3/n the ideal-gas compressibility 



(Bll) 



Finally, by inserting the solution of Eq.(B3) on the right-hand side of Eq.(3.4) and 
replacing Cj- ? (k, r) = (p//?)e _7T ^^ T with CV(k, r) in the mode-coupling integral, 
we immediately obtain Eq.(3.5). The derivation described here is entirely phe- 
nomenological. It can, however, be shown that the same result is obtained by using 
diagrammatic methods [1,22]. 
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